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Abstract 

Rare decay modes of the newly discovered standard-model-like Higgs boson h may test the 
flavor changing couplings in the leptoquark sector through the process h ^ pC. Motived by the 
recently reported excess in LHC data from the CMS detector, we found that a predicted branching 
fraction Br(/i —)■ pA) at the level of 1% is possible even though the coupling parameters are 
subjected to the stringent constraint from the null observation of r — py, where the destructive 
cancellation among amplitudes is achievable by fine tuning. 
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I. INTRODUCTION 


The newly discovered Higgs boson h at the mass 125 GeV is consistent with the Higgs 
boson predicted in the standard model (SM) [T]. The narrow decay width of a predicted size 
about 4 MeV in SM provides hope that the unusual rare decay due to new physics (NP) can 
have a measurable branching fraction. Recently, the CMS collaboration has reported [2] a 
possible excess in the decay process h —)■ with a signihcance of 2.4 a in the search for 

the lepton flavor violation (LFV). Assuming SM Higgs production, CMS obtained the best 
fit for the branching fraction summed over r“/i^ and , 

Br(h ^ = 0.84+°;p % . (1) 

We understand that it is too early to draw a positive inference until future analyses of 
higher statistics from both CMS and ATLAS experiments are performed. However, the 
present sensitivity at the 1% level is interesting enough to call for possible NP to deliver 
such a detectable rate but satisfy other rare decay constraints such as r —)■ /ly, 

Br(r —/ly) < 4.4 x 10“® at 90% C.L. from BaBar experiment [3] (2) 

at low energy. Indeed, there are a lot of theoretical activities [41122] along this line of inves¬ 
tigation. There were also a number of studies on LFV Higgs boson decays in literature [23] . 
We are particularly motivated by the leptoquark (LQ) associated with the third generation, 
which provides a large top quark mass insertion in the loop diagram. However, the LQ in¬ 
teractions also give rise to amplitudes for r —)■ py. We notice that the cancellation between 
two types of LQ contributions is possible for r —)■ py, leaving a large detectable decay rate 
for h —)■ r^/i^. Each of these two types of LQs has been outlined in the literature, such as 
in Ref.[l], but the combined version necessary for the cancellation was overlooked. 

The organization of the work is as follows. In the next section, we describe the LQ 
interactions associated with the top quark and tan lepton. In Sec. HI and IV, we calculate 
the decay r —)■ py and h —)■ respectively. We give details on numerical results in Sec. 

V. We conclude in Sec. VI. 
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FIG. 1. (a) Dipole transition of r —)• /i7 via the singlet LQ x. (b) Dipole transition of r —?■ ^^7 
via the doublet LQ D. 

II. LEPTOQUARK INTERACTIONS ASSOCIATED WITH THE TOP QUARK 

We associate the new LQs with the top quark of the third generation in order to avoid the 

very stringent constraints upon the flavor non-conservation among the first two generations. 

On the other hand, the mass insertion of the top quark can enhance the rate of the rare 

LFV Higgs decay mode among the second and the third generation leptons. To satisfy the 

1 

electroweak gauge symmetry, we can classify two types of LQs, the first one ys is a weak 

5 2 

SU{2) singlet, and the other one SU{2) doublet, i.e. iV = (Qs^Qs). The superscript 

denotes the electromagnetic charge number. LQs transform under the SU{3) color group as 

3 just like quarks. The relevant interactions are given by 

1 1 

9lx^ {Q3)hLr,L - gRX^tRTR 

+g''[Xf etRLr,L - g'^RQ3,LTRVt ( r <H- /r ) h.c. (3) 

The Feynman diagrams (for r —>■ /ly) that involve these LQ interactions are shown in Fig. 
The shown e symbol, basically ia^, links two SU{2) doublets into a gauge invariant singlet. 
For brevity, we do not show other Levi-Civita symbols that contract Weyl spinor indexes. 
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Also, (Qs)^ = (tL,bL), Lr^L = tl)'^ . The terms from exchanging t ^ fi are needed to 

induce the LFV between the muon and the tan. 


III. r ^ ^7 AMPLITUDES INDUCED BY LEPTOQUARKS 


We start with the contributions from the singlet x- We dehne to be the charge con¬ 
jugated state of t. In this way, we can avoid the use of the unfamiliar Feynman rule for 
two fermions flowing into a vertex. Instead, one fermion flows in and the other out. For 
example, the incoming r enters the hrst vertex and turns into a departing plus a boson 
X“5. The relevant vertices for the process r —)■ /ry are 

9lix~hH'^^L) - ^ ^ h ) + h.c. (4) 


For the outgoing left-handed muon, the Feynman amplitude that the external photon line 
attaches to the line is given by 


—)■ fi'j) = 


eqt-9ylrnt . 


(1 - zfdz 
Iq zmX + (1 — z)m: 




(5) 


where R stands for the right-handed chiral projection operator (1 -|- 7^)/2. It is understood 


that the external spinors ^(/i) and u{t) sandwich the Dirac chain. We keep track of the 
color factor 3 by a subscript c. Another amplitude where the photon attaches to is 

(1 — z)zdz 


, ^9^-19 R9^mt 

iM2{T -)■ /iy) = ^ -3n 




Iq zml + (1 — z)m: 




( 6 ) 


We set charges qt^ = —| and q _i = —Using z {1 — z) m. Adi, we obtain 


f Ad 1+2 — 


eg^g^ - |^(1 - 

IbTT^m^ Jq {1 — z)-\- zml/ 


The numerator of the integral becomes z‘^ — The overall result is 


fAdi +2 = (6(a;t) - |^o(a:t)) a^^KR , Xt = ml/m\ , and 


(7) 

( 8 ) 


^n{x) — 




Inx -|- (1 — x) -f ■ • ■ -|- 


(i-y) 


n+1 


n+1 


) l + {x-l)z 
and ^-i(a;) = 


dz 


(1 - xY+‘^ 

lux 


(9) 


1 -I- (x — l)z 


1 — X 
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So the amplitude is related to the integral function, 


Hi{x) =ii{x) - ^^o(a;) = [7 - 8x + + 2(2 + a;) ln(a;)] . (10) 

Our result is different from that in Ref.|l]. Note that there is another chiral amplitude for 
the outgoing right-handed muon, using These two amplitudes do not interfere in the 

zero muon mass limit. 

Now we switch to the contributions from the LQ doublet O. The relevant vertices for the 
process r —)■ /r 7 are 

-g''"^{m)(JlTR) + g'^L{^^)^^^L) + (rO/r)-h h.c. (11) 


For the outgoing left-handed muon. 


■ KA' ( ^ ^ eqtg'^g'^Lmt 

zAf i(r^/i7) =--3, 


(1 — zYdz 


'0 + (1 - z)m 




( 12 ) 


This corresponds to the diagram that the external photon line attaches to the t line. Another 


amplitude where the photon attaches to fl 3 is 


(1 — z)zdz 

Iq zrn^ + (1 — z)mi 




We set charges qt=\ and g^_5 = —|. Using z •H- (1 — z) in Adi, we obtain 




Ibvr^mQ 


A -§^2 _ - z)dz 

lo (1 - ^) + zmf/ml 




The numerator of the integral becomes — y. The overall result is 

IT Ifl 

iM\+2 = ^T* 3c (6(a^t) - 1 ^ 0 ( 2 :*)) cf^'^KR , x[ = ml/m 


IbTT^m^ 


So the amplitude is related to the integral function, 

1 


H2{,x) = ^i(a;) - 1^0(a;) = 


[-1 + 8a; - 7x^ - 2(2 - 5a;) ln(a;)] 


(13) 


(14) 


(15) 


(16) 


6(1 — a;)^ 

Note that there is another chiral amplitude for the outgoing right-handed muon, using g'/g'^ 
In general, the low energy effective operators of dim 5 are 


£efr D — [ha“^(C'z.T + CrR)t] + h.c. 

mt 


(17) 


Cr = 7^ (aR9L^tHi(x,) + g"Rg'lx[H2(x',)) 


( 18 ) 
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FIG. 2. h ^ Tfj, via the singlet LQ x- 


Cl = {glg^XtHi{xt) + g'lg'%^x[H2{x'^)) . (19) 

The partial decay width of the process r —)■ /ry is 

r(r /iy) = {\Cl\^ + \Cr\^) ■ (20) 

The general loop formulas for the radiative transitions can be found in Ref. [21] 


IV. h-fT + g VIA LEPTOQUARKS OF THE 3RD GENERATION 

For the rare decay h —)■ r/r, we start with the contribution from the SU{2) singlet 
_ 1 

leptoquark y 3 to the chiral amplitude of the outgoing right-handed r. We take the zero 
mass limit for g and r. At the one loop level, the Higgs coupling to t(pi)/2(p 2) is induced via 
a triangle diagram, which involves internal y lines. First, we concentrate at the diagram 
that the external Higgs touches the internal line, as shown in Fig. 


iM 


<1 _ 

x.R - 


3 , {-glgl) 


L{/+ ^2 + mt)L 

{2n)^ + - p2y - ml){P - m^) 


( 21 ) 


We use the Feynman parameterization trick to carry over the integration. The parameters 
a,/3,y are assigned to the denominator factors {i + pi)^ — m^, {i — P 2 Y — mf, C — rri^ 
respectively, under the constraint a + (3 + j = 1. Then we complete the square of the 
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denominator as follows, 


a[{t + pif - ml] + /?[(£ - p2f - ml] + - m^] 

= + 2api ■ t — 2/3p2 ■ ^ + ap\ + (3pl — {a + (3)ml — 7771 ^ 


= {i + api — l3p2Y — m?{a, f3) , where 777 ^( 0 ;, f3) = m^. — a/3s + (a + /3){ml — 777 ^) . (22) 

Shifting the loop momentum, we simplify the numerator of the the Dirac matrices with the 
equation of motion, 

(/^ + ml) —)■ — 2al3pi ■ p2 + ml —)■ — al3s + ml . 


Here the s variable is simply 2pi ■ p 2 = m\. The amplitude becomes 


7777 


iMln = - 3 = ( 9 « 97 “) 2Uad? j 


(23) 


The domain £\ covers positive a and (3, as well as a + /3 < 1. 

We perform the Wick’s rotation by Euclideanizing —)■ iqE 4 , —t —(^e^ 

d^QE —)■ d^q^dqEi = 47r|qr^pd|g^|dg£;4 —)■ 47r(g|; cos^ 0)|dg|;d0 —)■ n'^q'^dq'^. So 

Ml, = -33 («-) f 2'.dadff f . 

^ ' Ji, J -[qE+ m^{a,/3, s)f 16n^ 


9R9^mt 


A2 


3 aids —ml \ 

IbTT^n 7 ]^ V ^m‘^{a,/ 3 ,s) 2 2777^(0;,/d,s)/ 


(24) 


The logarithmic divergence has to be canceled by the one-particle reducible (IPR) diagrams 
with bubbles in the external lepton lines. The h line is either attached directly to r or p, 
picking up respectively the mass couplings mr or 777 ^, which are canceled by the propagators. 
It can be shown the corresponding IPR contribution to be 

A2 


+3/43^ t I log 


IbTT^n 


'0 


77772 -l- (1 — 7 ) 777 ' 


1 d'jL . 


(25) 


Overall, logA^ terms cancel. Therefore the combined amplitude is 

1 777 * 




-G, {9R9iL + 9%9lR) . 


(26) 
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Gx = 


1 1 al3s — mf \ 

m2(a,/?,s) “ 2 ^ 2m2(a,As) J 


log 


A2 


7m| + (1 — 7 )^? 


ci7. (27) 


Note that in the intermediate step, we can choose an arbitrary A for the convenience of the 
calculation. 

Alternatively, one can use the Passarino-Veltman [25] (PV) functions. The integral before 
the Feynman’s parameterization as given in (21) is 


I 




= I 


((^+Pi)^-m^) ( {l-p2)'^-m'f ) {P-rri ^) 


= / 


d^l 


m^+mf 


(27r)4 V (F+Pl)’^-”it)(F-P2)2-m?)p2_m2) + ((r+pi)2-m?)((r-p2)2-mf) 


(28) 


The hrst term simply gives the triangle function 


+ ml)Co{0,0, . 

The second term after shifting the loop momentum gives the bubble function 

-L_Bn(s m^) = f _^ r _ d-^l'I{2 -kY _ 

The result including the IPR bubbles is 


(29) 


(30) 


^ (ml + 0, s, ml 7nl 7nf) + i?o(s, mf, ) - i?o(0, ml ml . 


(31) 


We have cross-checked numerically that the G value from PV functions and from the Feyn¬ 
man parameterization method match each other. 

For the Feynman diagram that the Higgs touches the leptoquark, the required vertex 
originates from the bosonic interaction of The G coefficient is updated with 

the new addition, 

<Tx ^ G'x + C'o(0, 0, s, ml, m], ml) . 

When we come to the contribution from the SU{2) doublet leptoquark H, it is easy to see 
the simple translation, 

X3 77 H3 , 77 g'^!^ , 77 An , m^ ^ m^ , etc. 


Here m^ is the mass of the | charged leptoquark. More explicitly, 

= {ml+ml)Go{0,0,s,m1,ml,ml)+Bo{s,ml,m^)-Bo{0,m^,ml)+Xy^v^Go{0,0,s,ml,ml,ml) , 
















(32) 

rjl) = ~^cY^ — [{GxgRgL + Gng]^g"^)L + {G^g^gl + Gng'^g']^)R] . (33) 

The partial decay width, summing both processes h —)■ is 

r(h -)■ r'Fp^) = 20487r5^^ (^) + \GxgRgl + • (34) 

V. PHYSICS POSSIBILITIES 

To suppress the highly constrained r —)■ py, we tune the cancellation 

gWi^tHii^Xt) + g'^Rg'^Lx'tH2{x[) ^ 0 , (35) 

glgRXtHi{xt) + g'lg'%^x[H2{x[) ^ 0 . (36) 

We choose a simplihed scenario that only one chiral mode of the muon interactions is im¬ 
portant. Say, g^ 3> g^ and g'^ S> g'%, then we only hnely tune the corresponding one 
constraint, i.e. the first of the two. The ratio of the couplings gfigi/{g'fig'R} is given in 
Fig-i for the tuned cancellation in r —)■ /ly. The contour plot demonstrates that a large 
parameter space remains available for the required hne-tuning. 

Figure shows the predicted numerical size of Br(h —)■ both) versus g^i^gf for 

various LQ masses when the tuned cancellation is satished. We have set = 0 in our 
numerical study. A desirable branching fraction at 1% level occurs for the coupling product 
g]igf — 0.3 — 1 for the cases that tuq = from 600 GeV to 1 TeV. 

If we switch off either one of the canceling amplitudes in r —)■ /ry, the individual contri¬ 
bution to the Br(r —?■ /ry) is shown in Fig. This demonstrates how much hne-tuning is 
required. The Br(r —)■ /xy) would be at about the 1 % level for 500 GeV LQ and g^g^ about 
0.3 to 0.8 if only using one of the two canceling amplitudes. To go down from 10“^ to 10“® 
in the branching ratio, the two amplitudes are required to cancel each other by almost one 
part in 10^. 
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FIG. 3. Contour plot of the coupling ratio g^igi/{g'^Rg'i) on the plane, satisfying the 


tuned cancellation in Eq. (35) in the amplitude t ^ fij. 



FIG. 4. The predicted numerical size of Br(/i —?■ r^/i^ both) versus g^g^ for various LQ masses 
when the tuned cancellation is satisfied. The CMS Icr range of Eq. (|^ is also shown. 

Reference [1] also proposed a mechanism of cancellation in the amplitude of r —)■ /rq with 
the help of an additional vectorial top-like quark. However, the detailed gauge quantum 
numbers of the added structure have not been shown to be feasible. 

So far, we have set = 0. We show in Fig. [^the branching ratio Br(h —)■ for 
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FIG. 5. Individual contribution to the Br(r —)• 1 U 7 ) if either one of the two canceling amplitudes 
is switched off. 



FIG. 6 . The branching ratio Br(/i —)• both) versus for various choices of Xx ~ ~ 

—1, 0,1. The tuned cancellation is satisfied. 


various choices of = Xq 


-1,0,1. The tuned cancellation of Eq. (35) is satisfied. It 


gives additional freedom to achieve the desirable branching ratio for the rare Higgs decay. 
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VI. CONCLUDING REMARKS 


The rare decay of h —)■ can be at the current reachable sensitivity through the LFV 

LQ interactions, however hne-tuning is needed to avoid the stringent constraint from the 
null observation of r —)■ /ry. Here we have invoked more than one LQs, which couple to the 
third generation quarks, and the second and third generation leptons, in order to achieve a 
cancellation in r —>■ /ly but sizable contributions to h —)■ 

There is another issue related to possible contributions of these LQs to the muon anoma¬ 
lous magnetic moments (aka g — 2). It was shown a long time ago [26] and more recently 
[27] that if we choose, as we have chosen in the above analysis, the left-handed coupling to 
be much larger than the right-handed coupling for the muon, i.e. S> g^ and g'^ S> g''^, 
the LQ contribution to — 2 is highly suppressed by m^/MiQ and very small for the LQ 
mass range that we considered in this work. 

The required leptoquarks y and U can be strongly produced at the high energy and 
high luminosity hadron colliders. They have dominant decay channels into the top quark 
and the charged lepton r or /i. That is a very identihable signature. Both ATLAS [29] 
and CMS [281130] collaborations have searched for the third generation leptoquarks via pair 
production by strong interaction. The CMS have searched for the third generation LQ with 
electric charged —1/3 (similar to the of this work) decaying to a top quark and a tan 

lepton. They put a limit of 685 GeV at 95% CL on [28]. On the other hand, both 
ATLAS [29] and CMS [30] searched for the third generation LQ with electric charge —2/3 
decaying into a b antiquark and a tan lepton (similar to 0“^/^ in this work), and put a limit 
of 534 and 740 GeV, respectively. 

Therefore, there are still plenty of mass ranges for beyond 685 GeV and for 

beyond 740 GeV that one can directly search for a top or bottom quark with 
a tan or muon at the Run 2 of LHG-13. 
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